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A key problem in making precise perturbative QCD predictions is the uncertainty in determining 
the renormalization scale fi of the running coupling a 3 (/x 2 ). The purpose of the running coupling in 
any gauge theory is to sum all terms involving the /3 function; in fact, when the renormalization scale 
is set properly, all non-conformal p ^= terms in a perturbative expansion arising from renormaliza- 
^-^ tion are summed into the running coupling. The remaining terms in the perturbative series are then 

identical to that of a conformal theory; i.e., the corresponding theory with j3 = 0. The resulting 
scale-fixed predictions using the "principle of maximum conformality" (PMC) are independent of 
the choice of renormalization scheme - a key requirement of renormalization group invariance. The 

t ■ results avoid renormalon resummation and agree with QED scale-setting in the Abelian limit. The 

PMC is also the theoretical principle underlying the BLM procedure, commensurate scale relations 
between observables, and the scale-setting method used in lattice gauge theory. The number of 
active flavors nj in the QCD /3 function is also correctly determined. We discuss several methods 
for determining the PMC scale for QCD processes. We show that a single global PMC scale, valid 
at leading order, can be derived from basic properties of the perturbative QCD cross section. The 
D . elimination of the renormalization scale ambiguity and the scheme dependence using the PMC will 

rS—i not only increase the precision of QCD tests, but it will also increase the sensitivity of collider 

experiments to new physics beyond the Standard Model. 
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f*** " A key difficulty in making precise perturbative QCD predictions is the uncertainty in determining the renormaliza- 

tion scale fi of the running coupling a s ([i 2 ). It is common practice to simply guess a physical scale \i — Q of order 
of a typical momentum transfer Q in the process, and then vary the scale over a range Q/2 and 2Q. This procedure 
is clearly problematic since the resulting fixed-order pQCD prediction will depend on the choice of renormalization 
scheme; it can even predict negative QCD cross sections at next-to- leading-order [1]. 

The purpose of the running coupling in any gauge theory is to sum all terms involving the /3 function; in fact, 
when the renormalization scale /j, is set properly, all non-conformal /? 7^ terms in a perturbative expansion arising 
from renormalization are summed into the running coupling. The remaining terms in the perturbative series are 
then identical to that of a conformal theory; i.e., the theory with (3 = 0. The divergent "renormalon" series of order 
a™/3 n n! does not appear in the conformal series. Thus as in quantum electrodynamics, the renormalization scale fi is 
determined unambiguously by the "Principle of Maximal Conformality (PMC)" . This is also the principle underlying 
BLM scale setting [2] 

It should be recalled that there is no ambiguity in setting the renormalization scale in QED. In the standard Gell- 
Mann-Low scheme for QED, the renormalization scale is simply the virtuality of the virtual photon [3]. For example, 
in electron-muon elastic scattering, the renormalization scale is the virtuality of the exchanged photon, spacclikc 
momentum transfer squared /j 2 = q 2 = t. Thus 

w i-n(Mo) y> 



where 



n (M o) = nW ' n( y (2) 

v ' o; i-n(io) v ; 



sums all vacuum polarization contributions to the dressed photon propagator, both proper and improper. (Here 
II(i) = II(i, 0) is the sum of proper vacuum polarization insertions, subtracted at t = 0). Formally, one can choose 
any initial renormalization scale /ig = to, since the final result when summed to all orders will be independent 
of to- This is the invariance principle used to derive renormalization group results such as the Callan-Symanzik 
equations [4, 5]. However, the formal invariance of physical results under changes in to does not imply that there is no 
optimal scale. In fact, as seen in QED, the scale choice n 2 = q 2 , the photon virtuality, immediately sums all vacuum 
polarization contributions to all orders exactly in the conventional Gcll-Mann-Low scheme. With any other choice of 
scale, one will recover the same result, but only after summing an infinite number of vacuum polarization corrections. 

Thus, although the initial choice of renormalization scale to is arbitrary, the final scale t which sums the vacuum 
polarization corrections is unique and unambiguous. The resulting perturbative series is identical to the conformal 
series with zero /3-function. In the case of muonic atoms, the modified muon-nucleus Coulomb potential is precisely 
—Za(—q 2 )/q 2 ; i.e., /i 2 = — q 2 . Again, the renormalization scale is unique. 

One can employ other renormalization schemes in QED, such as the MS scheme, but the physical result will be 
the same once one allows for the relative displacement of the scales of each scheme. For example, one can compute 
the standard one-loop charged lepton pair vacuum polarization contribution to the photon propagator at photon 
virtuality q 2 using dimensional regularization. The result in MS scheme for spacelike argument q 2 — —Q 2 is 



f 2 /-l 

log -_MS_ = 6 dxx{l-x)\og "" e ' - 1^ — ^. (3) 
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At large Q 2 this is 



log^f=log%-5/3; (4) 

i.e., n 2 -^ = Q 2 e~ 5 / 3 . Thus if Q 2 >> 4m 2 , we can identify 

«Ms( e ~ 5/ V) = Oi G M-L(q 2 )- (5) 

The e~ 5 ' 3 displacement of renormalization scales between the MS and Gell-Mann-Low schemes is a result of the 
convention [6] which was chosen to define the minimal dimensional regularization scheme. One can use another 
definition of the renormalization scheme, but the final physical prediction cannot depend on the convention. This 
invariance under choice of scheme is a consequence of the transitivity property of the renormalization group [3, 7-9]. 
The same principle underlying renormalization scale-setting in QED must also hold in QCD since the n/ terms 
in the QCD j3 function have the same role as the lepton Nt vacuum polarization contributions in QED. QCD and 
QED share the same Yang-Mills Lagrangian. In fact, one can show [10] that QCD analytically continues as a 

function of Nc to Abelian theory when Nc — > at fixed a = Cpa s with Cp = -& — ■ For example, at lowest order 

@o = ~t~ (if^ 7 — l n /) — *■ ~3tt l n / a ^ ^ c = 0' Thus the same scale-setting procedure must be applicable to all 
renormalizable gauge theories. 

Thus there is a close correspondence between the QCD renormalization scale and that of the analogous QED process. 
For example, in the case of e + e~ annihilation to three jets, the PMC/BLM scale is set by the gluon jet virtuality, just 
as in the corresponding QED reaction. The specific argument of the running coupling depends on the renormalization 
scheme because of their intrinsic definitions; however, the actual numerical prediction is scheme-independent. 

The basic procedure for PMC/BLM scale setting is to shift the renormalization scale so that all terms involving 
the j3 function are absorbed into the running coupling. The remaining series is then identical with a conformal theory 
with j3 = 0. Thus, an important feature of the PMC is that its QCD predictions are independent of the choice of 
renormalization scheme. The PMC procedure also agrees with QED in the Nc — > limit. 

The determination of the PMC-scale for exclusive processes is often straightforward. For example, consider the 
process e + e~ — > cc — > ccg* — > cebb, where all the flavors and momenta of the final-state quarks are identified. The nj 
terms at NLO come from the quark loop in the gluon propagator. Thus the PMC scale for the differential cross section 
in the MS scheme is given simply by the MS scheme displacement of the gluon virtuality: [ip^c = e ~ (pt> +Ps) 2 - 

In practice, one can identify the PMC/BLM scale for QCD by varying the initial renormalization scale [Aq to identify 
all of the /3-dependent nonconformal contributions. At lowest order fi = j^ (ll/3JVo — 2/3n/). Thus at NLO one can 
simply use the dependence on the number of flavors n/ which arises from the quark loops associated with ultraviolet 
renormalization as a marker for fto- 

In QCD, the n/ terms also arise from the renormalization of the three-gluon and four-gluon vertices as well as from 
gluon wavefunction renormalization. 



It is often stated that the argument of the coupling in a rcnormalization scheme based on dimensional rcgularization 
has no physical meaning since the scale [i was originally introduced as a mass parameter in extended space-time 
dimensions. However, the QED example above shows that the MS scale is unambiguously related to invariants in 
physical 3+1 space. The connection of a-j^-g to the Gell-Mann-Low scheme can be established at all orders. This also 
provides the analytic extension [36] of the ctjjg scheme for finite fermion masses as well to timelike arguments where 
the coupling is complex. 

An example which shows how critical is to properly fix the renormalization scale is the thrcc-gluon vertex. The 
PMC/BLM scale which appears in the thrce-gluon vertex is a function of the virtuality of the three external gluons 
<7 2 ,<?f, and g|. It has been computed in detail in refs. [12]. The results are surprising when the virtualities are very 
different as in the subprocess gg — > g — > QQ. 

2 2 

fj 2 cx gmi ^ med (6) 
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where |g min | < |<z med | < |<7 ma xh i- e - 9max has the maximal virtuality [13]. The prediction based on simply guessing 
/i 2 ~ <lmax would give misleading results. 

The PMC/BLM scale that appears in the three-gluon vertex is the mass scale that controls the number of quark 
flavors nf which appears in the triangle graph. This is verified by keeping the quark masses and threshold dynamics 
in the loop. Thus we accurately determine the number of flavors nf that appears in the /3 function in the thrce-gluon 
coupling. This generalizes for all gluonic processes. 

Although these results have been obtained using the pinch-scheme, the final PMC/BLM result is scheme- 
independent. The pinch scheme is used because it provides a gauge-invariant setting for the analysis. In effect 
one calculates a scattering amplitude with three on-shell quark currents. One then obtains 14 invariant amplitudes 
which describe the three-gluon vertex, only one of which is renormalized. 

In fact the calculation of the PMC scale for the three-gluon vertex g a — > g\,g c given in Eq.(6) uses the pinch 
scheme to obtain a gauge invariant result. In effect, one computes the entire gauge invariant on-shell amplitude 
9a+<7a ~> <Z&<76 + <Ze9c including the triangle loop graph from quark loops with general mass. All 14 invariant amplitudes 
are computed analytically to one loop, only one of which is renormalized. The PMC scale for the three-gluon vertex 
as given in Eq.(6) also correctly sets the scale which controls the number of effective flavors which contribute to the 
/?— function for the three-gluon vertex. Details are given in refs [12], [13]. 

These results show that the usual method of guessing the renormalization scale for processes involving the three- 
gluon and four-gluon couplings, typically misses this essential physics, assigns rip incorrectly and mischaracterizes the 
perturbative prediction. The error which is introduced can be in principle eliminated at infinite order, but only if one 
can sum the renormalon series. 

The explicit result for the PMC/BLM scale is the physical scale controlling the quark threshold in the specific 
renormalization procedure used, but it is always possible to relate one scheme with another by the transitivity property 
of the rcnormalization group. This property is guaranteed by the PMC so there can be a constant displacement between 
schemes. 

The PMC method is a general approach to set the renormalization scale in QCD including purely gluonic processes. 
It is scheme independent and void of renormalon growth due to the absence of the /?— function terms in the perturbative 
expansion. We stress that the ^-function is gauge invariant in any correct renormalization scheme. The resulting 
conformal series is then gauge invariant. Thus the PMC is a gauge-invariant procedure. 

It is sometimes argued that it is advantageous not to fix the renormalization scale at all, since its variation provides 
a measure of higher-order contributions to the theory predictions. In fact, one obtains sensitivity only to the j3- 
dependent non-conformal terms by this procedure. In some cases the conformal contributions may be unexpectedly 
large. For example, the very large electron-loop light-by-light scattering contribution [14] ~ 18(a 3 /7r) 3 to the muon 
anomalous magnetic moment is unassociated with renormalization or the /3 function. Of course, one can still compute 
the variation of the prediction around the PMC scale as an indicator of higher order non-conformal terms. 

Stevenson has proposed that one should set the renormalization scale at a point where the predicted cross section 
has minimal variation with respect to /j, - the "principle of minimal sensitivity" (PMS) [15]. However, unlike the 
PMC, the application of the PMS to jet production gives unphysical results [16] since it sums physics into the running 
coupling not associated with renormalization. Worse, the PMS prediction depends on the choice of renormalization 
scheme, and it violates the transitivity property of the renormalization group [17]. Such heuristic scale-setting methods 
also give incorrect results when applied to Abelian QED. 

It should be emphasized that the factorization scale which enters predictions for QCD inclusive reactions is intro- 
duced to match nonperturbative and perturbative aspects of the parton distributions in hadrons; it is present even in 
conformal theory, and thus its determination is a completely separate issue from renormalization scale setting. 



II. IDENTIFYING THE RENORMALIZATION SCALE USING THE PRINCIPLE OF MAXIMUM 

CONFORMALITY 

Given the analytic form of the hard process amplitude or cross section as a series in a s (/Xg) evaluated at an initial 
scale /Uqi one can identify the PMC scale in a systematic way: 

1 . The variation of the cross section with respect to log fig can be used to distinguish the conformal terms versus 
the nonconformal terms proportional to the /? function. 

2. The identified nonconformal terms have the form j3 x log pij / '/i 2 , where p^ = pi ■ pj are the scalar product 
invariants i 7^ j which enter the hard subprocess. In practice, these terms can be identified as coefficients of rif, 
the number of flavors appearing in the /3 function; i.e., the flavor dependence arising from quark loops associated 
with coupling constant renormalization. The nj terms in QCD arise from the renormalization if the three-gluon 
and four-gluon vertices as well as from gluon wavefunction renormalization. 

3. The scale is then shifted p^ — > fi 2 in order to absorb the non-conformal terms. Thus when the scale is correctly 
set, the coefficients of a s (^ 2 ) become independent of the /3 function and log/x 2 . 

4. The series is then identical to that of the conformal theory where /3 = as given by the Banks-Zaks method [18]. 

5. The PMC scale is fixed for an observable (such as a differential cross section). PMC then can give a single 
effective global scale for the whole set of skeleton graphs entering the calculations which sums all the non- 
conformal /3-terms associated with renormalization into the running coupling. 

Other examples of this procedure will be given in the next sections. 

A. The Global PMC Scale 

Ideally, as in the BLM method, one should allow for separate scales for each skeleton graph; e.g., for electron-electron 
scattering, one takes a(t) and a(u) for the i-channel and u-channel amplitudes, respectively. 

Setting separate renormalization scales can be a challenging task for complicated processes in QCD where there are 
many final-state particles and thus many possible Lorentz scalars p?- = pi ■ pj. However, one can obtain a useful first 
approximation to the full PMC/BLM scale-setting procedure by using a single global scale /x 2 which appropriately 
weights the individual BLM scales. 

The global scale can be determined by varying the subprocess amplitude with respect to each invariant, thus 
determining the coefficients fy of logp 2 //x§ in the nonconformal terms in the amplitude. The global PMC scale is 
then 

^CxH;^, (7) 



i.e., 



log ii 2 = J2 w v lo S V% + log C (8) 



where the weight for each invariant is 



Jij 



2sijij Jij 



(9) 



and Y^i^-j w ij = 1- T ne constant C is the scheme displacement; e.g., C = e~ 5 / 3 for MS for /1 2 >> im*. 

As a specific example of the application of a PMC global scale, consider the electron-electron scattering amplitude 
in QED. (For simplicity, we will just take the contribution of the convection current to the amplitude, as in scalar 
QED.) The Lorentz invariant Born amplitude at the initial scale t is then 



M°f 



\t,u)=4na{to){^— + - ). (10) 

t u 

The running QED coupling a(q 2 ) in QED sums all proper and improper vacuum polarization graphs 

M(t,u) =4ira(t)(?— -) +47ra(w)(^— ) (11) 



where to leading order 

a (t) = a (i ) (1 + n,^l log ^). (12) 

V 37T to 

Aside from power-suppressed contributions involving the lepton masses, the resulting series is identical to the corre- 
sponding conformal theory with f3 = 0. 

In this process we have contributions from both the t - and u- channel amplitudes which require separate renor- 
malization scales for each skeleton graph. However, at leading order we can weight the amplitudes to obtain a single 
PMC/BLM scale which still sums the nonconformal /3 terms into the running coupling a(/i 2 ) at leading order. For 
example, using the standard Gcll-Mann-Low scheme, we can write 

M(t, u) = f(t)a(t) + g(u)a{u) = (f(t) + g{u))a(p?) (13) 

where f(t) = 47r(s — u)/t and g(u) = A-k(s — t)/u are the Born amplitudes for the t - and u -channels, respectively. 

Then in this case we have two basic PMC scales a(t) and a(u) for each skeleton graph in the standard Cell Mann- 
Low scheme used in QED. These couplings then sum all of the vacuum polarization corrections to the skeleton graphs 
to infinite order. The result is then gauge invariant and the logarithm of the global scale is 

log A 2 - .J®, r log(-*)+ .jy , r log(-u) (14) 

f(t) + g(u) f{t)+g(u) 

which duplicates the multi-scale result at NLO. 

One can also use the mean value theorem to obtain an effective single scale which analytically reproduces the exact 
multi-scale result to next to leading order. Since it matches the exact result at NLO, it also retains gauge invariance 
at this order. Moreover, the PMC single or multi-scale result is independent of the choice of scheme. The single scale 
result illustrates why it is wrong to guess a single scale like [i 2 = "pr since it fails to agree with this simple example. 

Using kinematical constraints such as the total momentum conservation s + t + u = the weighted scale dependence 
can be confined into the \og(t/u) term inside the running coupling. The global scale p, 2 is maximal at 9cm = t/2 
(/i 2 = \/tu = — t = —u) and vanishes at the boundaries (0,7r) where tan 2 (9cm /2) = t/u. The results are shown in 
Fig. 1. 

2t 




FIG. 1: The PMC/BLM scale as function of the CM angle 0cm :ee ->ee scalar QED 



III. A PMC EXAMPLE FOR QCD: APPLICATION TO JET CROSS SECTIONS IN 
ELECTRON-POSITRON ANNIHILATION 

As an example of the application of the PMC to QCD, we will show how the renormalization scale can be determined 

annihilation into two and three jets in AIS scheme. 



The two-jet cross section has only infrared divergences: 

~tf~> ^ "'"> r(2-A) 



(2) ,4tt/A vs/i Wo , T(l-A/2) 

^ =(J o(— H (i - A / 2 )-F7^ — rr ( 15 ) 



where a =47r^7V c Efiief. 



Here A = 4 — n is the number of extra space-time dimensions used to regulate infrared and ultraviolet divergent 
integrals. Eventually all of the infrared divergences and the factors involving A will cancel out. In dimensional 
regularization the scale \x is introduced as a mass scale to restore the correct dimension of the coupling. The gauge 
coupling g R is related to the renormalized coupling constant a R by 

9R _ a ^ ) (■ 2NA/2 7 gA/2 nfi N 

( 47r )(4-A)/2 - 47r V* J e ^ i0 ^ 

and here 7_g is the Euler constant. 

As discussed in the introduction, the mass scale of schemes defined by dimensional regularization attains its physical 
meaning when it is applied to QED. The renormalized gauge coupling is also related to the bare coupling by: 

9r = \f~Z^Z 2 /Z 1 g 07 (17) 

where Z\ is the renormalization constant for the quark-antiquark-gluon vertex, Z 2 for the quark field and Z3 for the 
gluon field. The renormalization constants are: 



*i = 1-^(^ + ^)1^—- — ) (18) 



A^n.(J—2-\ 



167T 2 \X UV XirJ 

5 ° ( 5 N 2 N\f 2 
f&F { 3 Nc ~ 3 Nf) ( X^y X IR 



Z* = l + q&(^»-^/)( — " — ) ( 2 °) 



where Ayy, A/# are related respectively to the UV— ultraviolet and IR— infrared poles. In the MS only the pole 
associated with UV renormalization is subtracted out, and this leads us to a redefinition of the gauge coupling: 

1 , 9% ,2„ 11 s 1 

g R 167r z 3 3 Xuv 

A suitable renormalization scheme is the A/5 which differs from MS by a constant term and the respective counterterm 
can be inserted in the Born cross section by shifting the coupling constant: 

^ = ^{l-(y^-^)^(i + (ln4 7 r- 7£ ))U«^{l- / 3 ar (1)1 (22) 



where: 



i = i + (ln47r- 7B ), (23) 



e e 



1 /ll 2 

A) = 7H [irN c - -T R (24) 



2tt V 6 3 

with T R = N f /2 , e = Xuv/2. 

The Born cross section for e + e~ — > q(j>i)q{p2)gipi) for massless quarks and gluons is 

da( 3 V A — 2 - 1 



(2) /47T/X 2 ,A/2 1 , >«f Mp n y-A/2S/ N 

= ^ U (^2-) r ( i-A/2) FA(xi '" 2) ^^ CFi3 (X1 '" 2) 



(25) 



Here 

Fx(x 1 ,x 2 ) = [(a* + a; 2 - 1)(1 - !!)(! - x 2 )]- A / 2 (26) 



and 



B v ~ x / 2S (x 1 ,x 2 ) = B v (x u x 2 ) - ^B s (x u x 2 ) 



(27) 



B v (x ll x 2 ) 



ry^ 1 I T*" 



(l-siXl-a*) 



(28) 



5 s (a;i,a;2) 



(l-Ki)(l-aj 2 ) 



(29) 



where Xi — — J= in the e + e CM. In terms of invariants: y.y = Sij/q = (pi + Pj) /q . Then x\ = 1 — y 2 3,x 2 = 

1 - J/13, »3 = 1 - j/i2, a?i + a; 2 + £3 = 2. 

The renormalized one-loop corrected cross section for e + e~ — > q{pi)q{p 2 )g{p3) is given by Eq. (2.11) of Fabricius 
et al. [21] For our purposes it is sufficient to quote only the term proportional to /3q in the M S— scheme: 



da& 



dx\dx 2 



da( 3 \n 2 ) 



Bioop dxidx 2 



. 2 .r(l-A/2),47T/Z 2 sA/2 _ ,. [i 2 , 

l + a s (fi 2 ) r(i_A) (^) ^( lo g^) + 



(30) 



where the coupling is defined as in Eq. 22: aMs(e los4,r ~ 7E /i 2 ) = a^jg(/i 2 ). The remaining contributions are indepen- 
dent of rif and /3o 

We can eliminate the non-conformal log-term proportional to /3 by shifting the rcnormalization scale ctMsin • ) m 
the Born cross section Eq. 25 

/ „2 

a s (^ 2 ) ~ a s (g 2 ) f 1 -a s (<j 2 )/3 log[— 

however, it is first convenient to shift the scale to // 2 — > (/j.q). 
Then 



da^ 



dx\dx 2 



d*W(u%) 



oncloop 



dxidx2 



Born 



, 2 . T(l - A/2) / 47r//§ v a/2 „ /, Mo\ 
l + a.0*S) r(l-A) ("^) ^o(log^) 



(31) 



Naively one could simply fix the scale to y/q 2 , but the 3-jet cross section will still be affected by IR divergences; in 
order to apply the PMC/BLM prescription we will first need to include the 4-jet contributions. 



IV. NUMERICAL SCALE FIXING 



The complete differential 3-jet cross section has been calculated by Fabricius et al. [21], and we quote here the 
results for the /3q— dependent terms: 



dV 3 )(e,<5) a s (q 2 ) 

(To — Of x 



dx\dx 2 



2vr 



(32) 



B v ( Xl ,x 2 ) 



. 1 - cos 8. , „ 2 13 
l-a s {q 2 )p (log( )+ logx 2 - — 



B s (x 1 ,x 2 )a s (q 2 )^-\+0(8 2 )) 



(33) 



where £3 = (2 — X\ — x 2 ) and 



da^ 3 \e,8) = da^+dcr^(e,8) 



(34) 



is the sum of the 3- and the 4-jet s contributions. The cancellation of the IR-poles is guaranteed by the KLN 
theorem [19, 20]. 

The variables (e, 8) arc small quantities introduced in the virtual amplitude in order to define the soft and collinear 
4-jet contributions to the 3-jet cross section. In particular these quantities refer respectively to the fraction of the 



total energy and to the cone opening angle which define the phase volume for a 3-jet event (for more details, see 
Ref. [21]). 

In order to extract the PMC/BLM scale we first work in the MS-scheme, fixing an arbitrary renormalization scale: 
jj 2 = //q. It turns out that /3q term of the 3-jet differential IR safe cross section has the form: 



d 2 a^(e,6) 
dx\dx 2 



er — ~ ^F X 

Zir 



(35) 



B v (x 1 ,x 2 ) 



, 1 — cos 6. „ , ,„ 13 , q 2 

+ 2 log (2- Xl -x 2 )- — + log \ 

2 3 n 2 



1 - a s (fj, ) (3 log( 



-B s (x 1: x 2 )a s ^ 2 )^+O(5 2 )) + ---. 

In principle we can extract information on the terms in this formula performing a detailed analysis of the dependence 
of the (3q— coefficient on the invariants. Performing a blindfold study we can single out the (3 — coefficient by means 
of the /3q— derivative of the whole cross section or either by the n/— derivative since: 



d/3 
Then we can factorize out the Born amplitude Eq.25: 



dj_ x dP± 
drif dnf 



d*W(ri) 



dxidx2 



Born 



d rfV 3 >(e,fr//g) 
dpo dx\dx 2 



-o»(A*o) (log( 

B S {X!,X 2 ) 

'2B v (xi,x 2 ] 



1 — cos S , 



•2 log (2 -a* -a*) 



13 .4.1 ? 



0(5 2 )) 



and at the first order approximation the PMC/BLM scale can be fixed numerically imposing: 



da( 3 V 



dx\dx 2 



d dV 3 )(e,<5;/i 2 ) 



dn 



f 



dx\dx 2 



nf-0 



(36) 



= 



(37) 



In the numerical procedure at NLO the analytic form of the cross section is not needed; one must only keep track of 
the appearance of number of flavors nj arising from loop diagrams involving renormalization. This procedure , which 
has been shown at NLO here, can also be iterated to higher orders in a s , by keeping track of the ny-terms entering 
the /3-function, leading us to an improvement of the accuracy of the PMC/BLM scale ^ 2 PMC - 

Following this procedure we can include all the non-conformal (3 terms into the running coupling constant for every 
physical process, setting the renormalization scale at the PMC/BLM scale without necessarily knowing the PMC/BLM 
analytic form. Thus we end up with a cross section which is formally equal to the corresponding conformal expansion 
with j3 = 0. In this particular case the PMC/BLM scale has the form: 



2 
HTMC 



q 2 (2 -xi -x 2 ) 2 —e 



2 B v (xi ,X2) 



(38) 



In this case the coefficient depends on the parton energies x\,x 2l on the angle parameter 5, and on the scale ratio q 2 / $ 
(all these quantities can be written in the form of Lorentz invariants). The different contributions to the coefficient 
can be also identified, term by term, by considering the most differential cross section (i.e. for the 3-jet case the triple 
differential cross section), by performing the derivative (or logarithmic derivative) with respect to the corresponding 
invariant, and then isolating the constant term. This procedure will be discussed in detail in the next section. 

V. THE PMC/BLM SCALE AS A FUNCTION OF THE JET MASS RESOLUTION PARAMETER 



As shown by Kramer and Lampe [16], one can define a QCD jet by defining a resolution parameter y ■ s as its 
maximal virtuality. The jet then consists of particles with total invariant mass squared smaller than y ■ s. Using this 
definition, we will perform the integration of the entire three-jet differential cross section, including real, da^' , and 
virtual, da^ s \ contributions in order to have a IR safe quantity. This gives a y-dependent integrated formula with 



/?o dependent terms which can be absorbed into the argument of the running coupling, according to the PMC/BLM 
prescription. 

The entire differential three-jet cross section [22]: 

_o^+ = / . dxT[l-x-z,x,z}a s (Q 2 )(l-p a s (Q 2 )(\og[x]+log[z}-?- )) 

00 ay J J 3 



a s (Q 2 )(T(y)-(3 a s (Q 2 )(C(y) + ....)) 



(39) 



T(y)a s (Q 2 ) (1 - /3 a s (Q 2 ) 21og[^-]) = T(y)a s (fi 2 BLM ); 



(40) 



where : a — era Gf Q 2 /2tt , s = Q 2 , x = y 13 , z — y 2 3, 



T[x 1 ,x 2 ,x 3 ] 



x l + x l + 2.ti(x 2 + a; 3 ) 



^2^3 



(41) 



and T(y), C(y) result from the partial integration of the LO- and NLO- terms of the 3-jet cross section (for more 
details see Ref. [22] [16]). 

Then in the 3-jet case, the BLM-PMC scale as function of the jct-virtuality y, has the analytic form: 



2 9 _5 i C(y) 

A = ^pmc/blm = s x e ^ rm 



(42) 



A plot of the PMC/BLM scale against y, the virtuality resolution of the jet, in e + e~ — > qqg is shown in Fig. 2. The 
result agrees with the BLM scale calculated by Kramer and Lampe in the MS scheme. The PMC/BLM prediction 
is scheme-independent; the specific value of the renormalization scale is rescaled according to the choice of scheme so 
that all results are commensurate. The PMC/BLM scale also accurately determines n/, the effective number of flavors 
in the /3-function. As is clear from the QED analog, the renormalization scale reflects the virtuality of the gluon jet; 
it thus must vanish when the resolution ys vanishes. As noted by Kramer and Lampe [16], the renormalization scales 
determined by the ad hoc PMS and FAC (Fastest Apparent Convergence) [23] procedures have the wrong physical 
behavior at y s — >• 0, since they become infinite /i 2 — > oo as the jet resolution and gluon virtuality vanish. 



0.5 r 
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0.00 




FIG. 2: The PMC/BLM scale, ixpmc (plane line) as a function of the jet resolution parameter y, for e + e 

For comparison, y/y is also shown (dashed line). 



qqg- 
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VI. PMC/BLM SCALE FIXING IN THE 3 JET CASE: THE COMPLETE DIFFERENTIAL CROSS 

SECTION 

In the case of the complete differential cross section; i.e., the most differential cross section for a given process 
without any constrained variables, the PMC/BLM scales depend on the number of flavors nf and on the independent 
invariants entering the process. In the case of the three jets, we notice that the cross section depends on the color 
and flavor parameters nf, Nc, Cf and on the kincmatical invariants S12, S13, S23 where the label 3 refers to the gluon 
momentum, and the indices 1,2 refer to the quark and anti-quark momenta. On the other hand, the nonconformal 
terms entering the running coupling depend only on the number of flavors n/ and on a reduced number of kinematical 
invariants. These terms can be identified by first varying the number of flavors nf and then the invariant Sij, whereas 
the constant term can be extracted by simply subtraction at the final step. Starting with the triple differential cross 
section for three jets, which is given by the sum of the singular part of 4-jet differential cross section da^ and the 
real 3-jet cross section da^ ( for more details sec Rcf.[22]): 



rf(T (s) + dcr (3) 
dz dy dx 



- a s (Q 2 ) 



5(1 - x - y - z){T[z,x, y] 



+ 24S!) CH .... ) + 54«!) A , o( .... ) 



2tt 



2tt 



a s (Q 2 )/3 log[x*y] 



a s (Q 2 ) 



^-F[z,y,x]} 



(43) 



with <7o = (To Cfs. For simplicity sake we are using the notation (z, x, y) for respectively the final state gluon-, quark-, 
antiquark-energy. In order to extract the first order terms related to the (3— function we can start performing an ab 
initio analysis of the cross section. We can first single out the /3o coefficient by means of the (3q— derivative, or either 
by the number of flavors nf— derivative, using Eq. 36 and then we can factorize out the Born amplitude: 



da^(Q 2 



dz dy dx 



1 



Bm 



n a s (Q 2 )d(3 



da^ + da® 

dz dy dx 



log [a: y] - - 



+ O(o s ), 



(44) 



da^(Q 2 ) 



dz dy dx 



Born 



a — ^— T[z,x,y\ 6(1 -x 



Finally, we can extract the weight for each invariant by taking the logarithmic derivative: 



d 



(do-W(Q 2 ) 



dlog(xi) \ dzdydx 



1 



_d_ ( da^+daW 

a s (Q 2 ) d/3o \ dzdydx 



(45) 



where x% = (x, y, z). The constant term can be identified by subtracting out all the logarithm terms from the /3o 
coefficient. Then at first order approximation in the coupling constant, the /ipA/p-scale for the 3-jct differential cross 
section has the analytic form: 



H 2 PMC ~Q 2 xCx ]\x^=Q 2 



x y e 3. 



(46) 



A. Commensurate Scale Relations 



Relations between obscrvablcs must be independent of the choice of scale and renormalization scheme. Such 
relations, called "Commensurate Scale Relations" (CSR) [24-26] are thus fundamental tests of theory, devoid of 
theoretical conventions. One can compute each observable in any convenient renormalization scheme, such as the MS 
scheme using dimensional regularization. However, the relation between the observables cannot depend on this choice 
- this is the transitivity property of the renormalization group [3, 7-9]. For example, the PMC relates the effective 
charge a 9l (Q 2 ), determined by measurements of the Bjorken sum rule, to the effective charge a/{(s), measured in the 
total e + e~ annihilation cross section: [1 — a 9l (Q 2 )/tt] x [1 + an(s*)/ir] = 1. The ratio of PMC scales y/s /Q ~ 0.52 
is set by physics; it guarantees that each observable goes through each quark flavor threshold simultaneously as Q 2 
and s are raised. Because all j3 7^ nonconformal terms are absorbed into the running couplings using PMC, one 
recovers the conformal prediction [25]; in this case, it is the Crewther relation [27-31]. Thus by applying the PMC, 
the conformal commensurate scale relations between observables, such as the Crewther relation, become valid for 
non-conformal QCD at leading twist. 
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VII. CONCLUSIONS 

As we have shown, the principle of maximal conformality (PMC) provides a consistent method for setting the 
optimal renormalization scale in pQCD. The PMC scale is determined by identifying the ft terms in the next-to- 
leading contributions and making the appropriate shift in order to include the /3-tcrms into the running coupling. 
This can be done most simply by identifying the nj terms which come from quark loops of skeleton graphs. This 
includes the nf terms which renormalize the three and four gluon couplings. This procedure has been used to identify 
the correct PMC scale for the threc-gluon vertex [12] [13]. The resulting series is identical to that of the corresponding 
conformal theory with ft = as given, for example, by the Banks-Zaks method [18]. 

The global PMC renormalization scale is particularly useful for very complex processes; one only requires the 
dependence of the calculated subprocess amplitudes on the initial renormalization scale /Xg and n/, the number of 
quark flavors appearing from quark loops associated with renormalization. The single global PMC scale, valid at 
leading order, can thus be derived from basic properties of the perturbative QCD cross section. 

We have discussed specific methods for efficiently determining the PMC renormalization scale analytically or numer- 
ically for QCD hard subprocesses. The analytic form of the PMC renormalization scale can be determined by varying 
the subprocess amplitude with respect to each invariant, thus determining the coefficients fy of logp|//Xo m the 
nonconformal terms in the amplitude. This result can be used to fix the renormalization scales for each contributing 
skeleton graph. However, we have shown that a single PMC global scale can then determined at NLO by appropriate 
weighting. Alternatively the numerical value of the PMC scale can be determined without specific information on the 
analytic form from the n/-derivative of the cross section. The two methods give rise to the same results at NLO. 

The factorization scale, in contrast, is the scale entering the structure and fragmentation functions. Unlike the 
renormalization scale, a factorization scale ambiguity occurs even in a conformal theory. The factorization scale 
should be chosen to match the nonperturbative bound state dynamics with perturbative DGLAP evolution. This 
could be done explicitly using nonperturbative models such as AdS/QCD and light-front holography where the light- 
front wavefunctions of the hadrons are known. 

Note that one applies the PMC method to renormalizable hard subprocesses (including the associated radiation 
diagrams required for IR finitencss) which enter the pQCD leading-twist factorization procedure. The initial and 
final quark and gluon lines are taken to be on-shell so that the calculation of the hard subprocess amplitude is gauge 
invariant. Thus the application of the PMC to hard subprocesses docs not involve the factorization scale, and thus 
no double or single logarithms which involve the factorization scale enter. 

The usual heuristic method of guessing the renormalization scale and varying it over a range of a factor of two gives 
scheme-dependent results, leaves the non- convergent perturbative series and gives the wrong result when applied to 
QED processes. In fact, varying the renormalization scale around such a guess only exposes nonconformal contributions 
involving the ft function; it gives no information on the conformal contributions. The PMS method [15] has similar 
faults - it violates the transitivity property of the renormalization group, depends on the choice of scheme, is wrong 
for QED, and as shown by Kramer and Lampe [16], leads to unphysical results. In contrast, the PMC method, 
which has no such disadvantages, and satisfies all principles of renormalization theory, gives the optimal prediction 
for pQCD at each finite order. 

The PMC is the theoretical principle underlying the BLM procedure and commensurate scale relations between 
observables - the rigorous scale-fixed scheme-independent relations in QCD between observables, such as the Gener- 
alized Crcwther relation; it is also the scale-setting method used for precision determinations of a s in lattice gauge 
theory [32]. In addition, it has been recently shown that for certain observables in 2-jet production the results of the 
MOM-BLM method arc very similar to those of MSYM theory [33] [34] [35]. 

In the case of the BLM method, one deals with separate renormalization scales for each skeleton diagram, as is 
done in QED. The PMC method provides a single effective renormalization scale which reproduces the BLM scales at 
NLO, even for rather complex processes that are in our list of important projects, such as VF+Jets, e + e _ annihilation, 
it production, and for general observables; e.g. differential cross sections, asymmetries. 

If one considers a process with high multiplicity, then one confronts separate BLM scale for each of the multiple 
skeleton diagrams; thus the number of BLM scales will appear as the jet multiplicity increases. The PMC method 
replaces these multiple scales with an effective single scale at NLO. 

We have discussed in this paper an illustration of the PMC procedure for 3-jet production in e + e~ annihilation 
where the nf terms arise from the inclusive 4 jet cross section after IR cancellation: these terms are included in the 
PMC scale with the effect of lowering its value. 

The PMC method provides the correct renormalization scale from first principles without ambiguity or renormal- 
ization scheme dependence. The residual errors from the resulting conformal series provide an accurate assessment 
of higher order errors. The PMC/BLM uncertainty is zero at the order computed. The PMC is equivalent to the 
standard method used to eliminate the renormalization scale ambiguity in precision tests of QED. 

The PMC method gives results which are renormalization scheme independent at each finite order. The PMC also 
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determines the correct number of flavors rij ; this is particularly important when one uses a renormalization scheme 
which is analytic in the quark masses such as the analytic extension of the MS scheme [36]; one can then include 
the correct flavor threshold dependences and transitions as one evolves the QCD coupling. The correct displacement 
between the argument of the schemes is also automatically determined. 

We stress that PMC does not capture all higher-order effects. One still has higher order corrections in the conformal 
series. These can never be discovered by varying the renormalization scale, since this variation only exposes terms 
proportional to the /3— function. It is incorrect to require the scale choice to remove all higher order terms. For 
example, in QED, the muon anomalous moment receives a large contribution at order a 3 from the electron-loop 
light-by-light insertion. This is due to the physics of the higher-order processes — not the running QED coupling. It 
is thus incorrect to vary the renormalization scale to minimize the effect of higher order corrections, since the variation 
of fin cannot expose large terms in the conformal scries. Thus the PMC correctly and unambiguously exposes higher 
order terms which are intrinsic to physical effects, unrelated to the QCD running coupling. 

We emphasize that the PMC method for setting the renormalization scale gives predictions for observables which 
are independent of the choice of renormalization scheme - a key requirement for a valid prediction for a physical 
quantity. The argument of the running coupling in a given scheme which appears in the resulting conformal scries 
has the correct displacement so that the result is scheme-independent. The number of active flavors n/ in the QCD (3 
function is also correctly determined, and the renormalization agrees with QED scale-setting in the Nc — > Abelian 
limit. Furthermore, the resulting conformal series avoids the need for renormalon resummation. 

A consistent application of the BLM/PMC procedure to B-decays, including B — >• X s + 7, has been developed 
including resummation to all orders in the strong coupling constant. A review and extension of this procedure is given 
by Melnikov and Mitov [37] 

The PMC procedure has recently been extended to the four-loop level, [38] demonstrating that it provides a con- 
sistent, systematic and scheme-independent procedure for setting the renormalization scales up to NNLO. 
The explicit application for determining the renormalization scale of R e + e - (Q) up to four loops has also been pre- 
sented [38]. The PMC procedure has recently been applied to the ti hadroproduction cross section, a major test of 
the standard model [39, 40]. The PMC scales are found to be smaller than used in the conventional scale-setting 
method. As a result, a larger a t i is obtained in comparison to the conventional scale-setting method, which agrees 
well with the present Tevatron and LHC data. 

Clearly, the elimination of the renormalization scheme ambiguity using the PMC will greatly increase the precision 
of QCD tests and increase the sensitivity of measurements at the LHC and Tevatron to new physics beyond the 
Standard Model. 
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